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Abstract. Solids with spatially varying photonic structure offer gaps to light of a wider range of frequencies 
than do simple photonic systems. We solve numerically the field distribution in a cholesteric with a linearly 
varying inverse pitch (helical wavevector) using equations we derive for the general case. The simple idea 
that the position where the Bragg condition is locally satisfied is where reflection takes place is only true 
in part. Here, reflection is due to a region where the waves are forced to become evanescent, and the 
rate of variation of structure determines over which distance the waves decay and therefore how complete 
reflection is. The approximate local Bragg-de Vries schemes are shown to break down in detail at the edges 
of the gap, and an analytical estimate is given for the transmission coefficient. 
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1 Introduction 

Photonic band structures can be made to vary spatially by 
varying the lattice repeat distance as one penetrates the 
sample. An example where this technique has been used 
with great success is that of cholesterics. In such liquid 
crystals the director (the direction of high refractive index) 
rotates on moving perpendicularly to itself to give a helical 
structure, see the sketch fig-H These structures are known 
to have a single photonic stop gap for normal incidence 
of circularly polarized light of the same handedness as the 
helical cholesteric and no gap for light of the opposite 
handedness. 

We study here variable pitch cholesterics (VPCs). The 
motivation for spatially varying structures is to have re- 
flection over a broader range of frequencies than simply 
that provided by the gap of a uniform system. This has 
been usefully employed in wide band reflectors and circu- 
lar polarizers for display purposes. The idea employed to 
interpret such effects is that light of a given frequency pen- 
etrates until the repeat distance of the solid is locally that 
of the wavelength of this light in the solid. At this point 
a Bragg-de Vries condition is satisfied and refiection oc- 
curs. Our purpose is to explore whether such an approach 
has any utility and what limitations arise, for instance 
unwittingly allowing unwanted transmission under some 
circumstances. The counter argument to this simple argu- 
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ment would be that Bragg effects are coherent - one needs 
very many such repeats to get sharp and total Bragg inter- 
ference and in a varying medium this is achieved nowhere. 
To some extent the answer must be that significant evanes- 
cence of the wave must occur between first encountering 
the gap and leaving it. We explore this numerically and 
by approximate methods. A rather different theoretical 
approach has used the Berreman matrix formulation 

Since cholesterics are locally nematic and therefore also 
fluid, a VPC must be solidifled to stabilize the pitch gradi- 
ent once it has been set up. One method is to photo- induce 
a concentration gradient in a blend of chiral and achiral 
photo-crosslinkable mesogens A network (elastomer) 
is formed which freezes in the structure. Otherwise, ther- 
mal diffusion between two cholesteric layers with different 
pitches can be used to set up a variation in pitch which 
can then be frozen in by quenching to the glassy state |S] . 



2 Governing equations 

In conventional cholesteric liquid crystals, the nematic di- 
rector n rotates uniformly in a plane as one advances in a 
line perpendicular to this plane. We shall adopt the follow- 
ing convention: the nematic director n is conflncd to the 
ccy-plane, and the direction in which the director changes 
is therefore the z-axis, see flg. ^ The director of ordinary 
cholesteric liquid crystals subsequently becomes 

2tt 

n = (cos (/)(z), sin(/)(z), 0) with (t){z) = — z, (1) 
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Fig. 1. A classical, constant pitch cholesteric structure. The 
arrows indicate the nematic director prevalent in the xy-p\ane 
at a certain z coordinate. Note that n and —n describe the 
same physical state, namely an orientational order along the 
direction of n (or — n). 



where pq is the pitch length. 

To allow for the most general phase evolution, we look 
at the governing equation for electromagnetic waves in a 
cholesteric liquid crystal whose pitch does not need to be 
a constant. In other words, we do not assume the specific 
form — (2t: /pq)z in eq. (^3), but leave the phase in the 
most general form (/'(z). On the other hand, we restrict 
ourselves to electromagnetic waves which propagate along 
the z-axis only, i.e. perpendicular to the director n. There- 
fore the electric and magnetic fields of the wave only have 
X and y components 

The anisotropy of the liquid crystalline phase described 
by the director n implies that locally the dielectric tensor 
respects the uniaxial symmetry of the liquid crystalline 
phase. It therefore assumes the following form: 



e ei^ + (£2 - ei) nn'"', 



(2) 



where ^ is the 3x3 unit matrix and where €2 and ei are 
the dielectric constants parallel and perpendicular to the 
director n respectively. Since n rotates in the xy-plane 
along the z-axis eq. the dielectric tensor e eq. |(2J) 
depends too on the z-coordinate. 

In the present case, the electromagnetic wave equation 
(assuming a unit magnetic permeability = 1), reads [Q: 



1 d^Dt 

c2 dt^ 



(3) 



with the dielectric displacement field Dt = eEt . 

Introducing the time dependence Et{z) = e*"*£^(z) 
(and similarly for D), all time derivatives become a mere 
multiplication by icu, and we obtain: 



d^E 



with ko 



(4) 



Note that the tensor e is in general not diagonal in the 
fixed laboratory frame (xyz) (cf. eqs. ^ and |2Jl), however 
due to its symmetry given by eq. the dielectric tensor 
is diagonal in the frame defined by the rotating director 
n(z). In this frame, which we denote by a tilde and (^,?y) 
coordinates, the electric displacement field is simply given 
by: 



We take rj to be the direction along n and hence typically 
£2 > ei- 

The coordinate transform into the rotating frame (at- 
tached to the director n), reads as follows: 



cos(/)(z) sin(/)(z) 
- sin (/'(z) cos (j){z) 



E = Eiz)E with E{z) = 
In this new frame, eq. (@J transforms into 



k^eE 



-dl + <i>' 2q^'d, + 0" 



(6) 



(7) 



where dz = d/dz and (p' — dip/dz, and e is now the di- 



agonal 2x2 dielectric tensor ^ J in the (^,77) frame. 
Dropping the tilde, here and henceforth, eq. lO reads: 

kl ei = -4' -I- E^(t)'^ + 2E'^4)' -(- E^4)" 

kl £2 Er, = -E',; + E,,4>'^ - 240' - E^4>". (8) 

3 Cholesterics with constant pitch 

Ordinary cholesteric liquid crystals are characterised by a 
constant pitch po- Hence the director n performs a steady 
rotation as one travels in z-direction. The physical system 
is periodic in although eq. is only periodic in pq, 
reflecting the fact that both n and —n describe the axis 
of anisotropy of a uniaxial system equally well, see eq. 0. 

For constant pitch, the phase 0(z) is given hy (p — qqz 
with 90 = 1^ = const, and hence (p' = qo and (f>" = - 
the classical problem of de Vries. The wave equation © is 
exactly solvable 1 in the (^, 77) frame using the solution: 



E 



A 
iB 



(9) 



We obtain two linear homogeneous equations for the coef- 
ficients A and B. There can only be a non-trivial solution if 
the determinant of the corresponding matrix equals zero: 

k^-e{kl[ei + e2] + 2ql) + [e^kl^ql){e2kl-ql) = . (10) 

This quadratic equation in k^ establishes a relation- 
ship between the temporal frequency lu — koc and the 
corresponding wave vector k. In the reduced zone scheme, 
the relationship fcg as a function of k is shown in fig.[5fa). 

The ratio A/B = E^/E^i of the field components also 
varies as fc(fco) is varied: 



2kqn 



9o - ^2kl 



(11) 



eiE^ and 



£2E„ 



(5) 



where k and kg are related by eq. (|10|l , see fig. IHb) . The 
choice of the dielectric constants ei and £2 in this illustra- 
tion corresponds to refractive indices perpendicular and 
parallel to n of ni = 1.5 and 712 — 1.8. 
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Fig. 2. (a) Reduced temporal frequency ko — cj/c vs. wave vector k: in a cholesteric liquid crystal with constant pitch po, 
there are two pairs of solutions. One of them exhibits a single stop gap for normal incidence at frequencies ko — uu/c between 
[9o/\/e2i ^o/x/eT] (solid line) with qo being the cholesteric's wave vector 2n/po. The other is gapless (dashed line). The plot can 
be extended to negative values of k by reflection in the feo-axis. In this example, we arbitrarily chose ei — 2.25, £2 ~ 3.24 (the 
dielectric constants perpendicular and parallel to the local director) and qo — 1. (b) The corresponding ratio of the flelds along 
and perpendicular to the rotating director as a function of frequency ko = cu/c. The dotted line of the inset shows how the ratio 
A/B becomes imaginary in the gap and that the variation of A/B is strong near the edges of the gap. 



Equation (|10|l allows two pairs of solutions: for each 
fixed frequency kg = ^, there are two solutions whose 
wave vector k have the same magnitude, but opposite sign, 
reflecting the fact that there are waves travelling in the 
two opposite directions along the z-axis. Apart from this 
trivial pairing, there are two genuinely different solutions 
to eq. H10|l . One of them exhibits a band gap in the fre- 
quency range uj/c G [<}o/ Qo/ Note that due to 
the periodicity in the cholesteric structure, the graphs in 
fig. Ufa) were shifted by 1 along the fc-axis to yield the 
reduced zone representation. 

The non-trivial pair of solutions, the de Vries eigen- 
states, are in general of elliptical polarization, the prin- 
cipal directions of which rotate with n. At the two band 
edges at ko = and qo/y/ei, the states are plane 

polarized with the plane of polarization rotating with n. 

Light with the same handedness as the cholesteric has 
states {E^,0) and (0, Ejj) at the band edges, corresponding 
to E respectively along n or perpendicular to n. The ratio 
A/B is then 00 or respectively. The former is along a 
polarizable direction and hence produces a lower energy 
that the latter even though the spatial variation k is the 
same — a gap has been produced. Light rotating with 
the opposite handedness does not couple coherently to the 
cholesteric and retains a gapless dispersion relation, the 
dashed line in fig. El The de Vries problem, in this photonic 
band structural language of fig. 13 can be extended to non- 
uniform but periodical director rotation ,7,. 



4 Cholesterics with varying pitch 

In varying pitch cholesterics (VPCs) the nematic director 
n rotates non-uniformly along the z-axis, i.e., the phase 
4>{z) of eq. 1^ does not have the simple form (f>{z) = qqz — 
2'Kz/po with (70 — const, but assumes a more general form. 
It has been reported that the local pitch can vary expo- 
nentially with z 12], i.e., we can write q — go exp(— kz). 



Other careful, scanning electron microscopy investigations 
showed the pitch to be linear over a range of z- 
values. Here, we assume the simplest possible inverse pitch 
variation, namely a linear dependency on z, leading to a 
quadratic dependence of the phase z: 

q = qQ+iqlz =^ (j){z) = {qQ + \-iqlz)z . (12) 



4.1 Local de Vries approximation 

A simple framework for considering a wide band circular 
polarizer, filter or refiector now presents itself As 
light of a given frequency enters and accordingly has a 
given wavelength in the medium, then since q{z) changes, 
there will be a depth at which a Bragg-like condition will 
be satisfied, and at this depth the wave will be totally 
reflected. This is where the wavelength in the medium 
matches the local periodicity 27r/q(z). We examine how 
much credence to put in this description. For Bragg re- 
flection to occur, one requires a semi infinite medium, at 
least a thickness of sample corresponding to very many 
wavelengths. In the variable pitch cholesteric, the Bragg 
condition is met strictly only at a point and a Bragg mech- 
anism fails. Indeed, instead of a total reflection, there may 
be circumstances where it is only partial, as we flnd below. 

One can take an analogous approach of applying de 
Vries theory locally to argue for reflection of light of a 
given wavelength at the particular point in the VPC where 
for light at this frequency a de Vries gap is encountered. 

In the constant pitch case, one can determine an eigen 
wave vector k given a frequency lo (or equivalently fco), 
see flg. 121 The pitch wave vector go played the role of a 
parameter in the function k — k(uj). Let us now change 
perspective and consider the frequency w or fco = ^ as a 
parameter, but allow the inverse pitch length q to change 
according to eq. (|12|l . In the following, we choose ko = 
!£ = 2,ei = 2.25, £2 = 3.24, qo = 1, and 7 ^ 0.05. 
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Fig. 3. Two of the four branches of the local wave vector fc vs. 
2-coordinate: assuming a fixed frequency fco, eq. links the 
wave vector k and the inverse pitch length go- Previously, go 
was taken to be a constant, but eq. 11211 introduces a depen- 
dence on the z-coordinate. Thereby, the relation between the 
wave vector k and the local inverse pitch length q becomes a 
function of z. For this plot, we have chosen ei = 2.25, £2 = 3.24, 
fco = 2, go = 1, and 7 = 0.05 (go7 being the spatial rate of 
change of cholesteric wavevector) . With this choice of param- 
eters, the frequency fco lies below the forbidden frequency range 
at 2; = 0. The forbidden zone of branch (2) extends over the 
region z G [zi, 22] = [40, 52]. 

By this reasoning, we are able to determine a wave 
vector fc satisfying the de Vries condition eq. (|10|l given 
a particular cholesteric inverse pitch length q, which in 
turn depends on the z-coordinate by means of eq. p2(l . 
Hence, we can plot an hypothetical wave vector k against 
the local position z, fig. El Note that for a given frequency 
A;o, there are again four solutions derived from eq. (|10l) . 
Figure 121 only shows two branches, one of them becoming 
imaginary between z\ and Zi, with Zi — {^/eiko — go)- 

The second pair of branches can be constructed trivially 
by introducing minus signs: k 1— > —k. In the region [zi, Z2] 
where the frequency fcg falls into the band gap for two of 
the four branches, the corresponding wave vector k be- 
comes imaginary, indicating an evanescent behaviour of 
the wave. We expect there a decaying and increasing so- 
lution (see eq. lO for imaginary k). 

One can now return to the question of reflection, or 
equivalently transmission. The band gap starts at zi and 
extends to Z2- Transmission depends on how close Z2 is to 
zi and how rapid the exponential decay in this gap is. 



4.2 Numerical investigation 

The above tentative interpretation of the solution being 
locally de Vries-like is qualitatively confirmed by looking 
at a full numerical solution in fig. 0] Here, we have cho- 
sen the initial conditions at z = such that we obtain a 
decaying solution in the forbidden zone as opposed to an 
increasing solution, which one obtains for different initial 
conditions. From the simultaneous field equations for 
and i?,,, it is clear that one could divide through by a 
constant, say £'j,(z = 0). The required 4 boundary condi- 




Fig. 4. Numerical solution of eq. Q with 4>{z) = {qo + ^'yqQz)z: 
the initial conditions a,t z — are chosen such that the solution 
decays over the evanescent range z G [21,22] = [40,52]. The 
numerical parameters are the same as in fig. El and the initial 
conditions as follows: _E^(0) = -0.036, £^(0) = -0.58, E'^{0) = 
2.0, and -E^(O) — —0.067 (the fields and their gradients in the 
rotating frame (f , 77) at 2 = 0). 

tions then reduce to 3 choices, the 4th being thus set to 
-E^ = 1 at z = 0, reflecting the arbitrariness of the overall 
scale in the problem. 

Several features of this solution are important for what 
follows. The amplitudes on either side of the gap region 
differ considerably. The small amplitude to the right of the 
gap after partial transmission gives a small Poynting flux, 
E AH, of energy in the positive z direction. Although the 
amplitudes on the left are much larger, they correspond 
to the sum of an incident and reflected wave resulting in a 
correspondingly small Poynting vector. The ratio A/B of 
the field amplitudes is approximately constant away from 
the band edges at zi and Z2. In the de Vries, constant 
pitch solution this ratio is respectively 00 and at these 
points, extreme values that are not achieved in this case 
of partial barrier penetration. The ratio also becomes very 
large at z = —20 which for our value 0.05 of 7 is the point 
where there is no rotation of the director at all, that is 
where g(z) — > 0. Locally the medium is achiral and there 
is no need for the two-component, de Vries-like solution. 

One can now qualitatively analyze reflection and trans- 
mission. With rising q vector on penetrating the VPC, one 
must eventually reach what would be the lower band edge 
in a uniform de Vries system. However, the gap translates 
into a spatial region of decaying gap states, and ultimate 
transmission is determined by how wide this region is. In 
particular does significant decay of the wave function oc- 
cur before the propagating region of z > Z2 is attained? If 
not, then transmission is higher and refiection lower. De- 
cay depends therefore in part on the gap's spatial extent 

Z2 - zi = ( - V^i ) (13) 

over which decay occurs. Clearly, if the spatial variation 
is weaker, i.e. 7 small, then the gap is spatially more ex- 
tended and decay has a greater length over which to take 
place (see also the discussion of the transmission coeffi- 
cient leading up to eq. 1)19(1 ^. 
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4.3 Approximate solutions 

For constant (inverse) pitch, the solution eq. © was a 
plane wave. In the case of a cholesteric with a z-dependent 
(inverse) pitch, we choose a more general trial function, 
namely a WKB-like (Wentzel, Kramers and Brillouin) ansatz 
assuming for simplicity that the two field components 
and E^i share the same phase S{z): 

Inserting eq. ((njl into the wave equation ((Sj) yields two 
linear equations for A and B, which have a non-trivial 
solution only if the corresponding determinant is zero. 
Therefore, the analogous equation to eq. I|10() is: 

+ S"j\l (z(ei + e2)fco' - - QoS") 

+S'^h' + S'A^qW (15) 
-S'^^'ql ((ei + e2)kl + 2qlw^ + 2iql^^S") 

£162^0 - (ei + (■2)klqlw'^ + q^w'^ + (?o7^ = 

where we have used a rescaled variable w = 1 + qo^z. 
Also note that the prime denotes the differentiation with 
respect to the new variable w: S' = dS{w)/dw. 

This equation is not just an algebraic equation like 
eq. (Il()|l . but rather a non- linear second order differential 
equation for S. An approximation is buried in this proce- 
dure. The components A and B were taken as constants, 
any variation with z they might have in common being 
taken up as logarithmic terms in the general phase S{z) 
(including any usual pre-exponential terms like (S")^^/^ 
that one would find in conventional WKB). However, we 
know from eq. (|ll|l and fig. El^b) that the ratio A/B in 
the uniform case depends on k which is now varying with 
respect to q{z) since the pitch is spatially varying. Hence 
even in a local de Vries picture, A/B must be a function of 
z through q{z). This variation has been ignored on insert- 
ing (|14() into eq. ((HJ. The full numerical solution, fig. 01 
indeed shows the ratio A/B to vary spatially. From fig. 
|2Ib), one can see that the de Vries value is oo at the lower 
edge and at the upper. However, the variation A/B goes 
to zero as 7 ^ 0, as expected, and can be estimated for 
small 7 to be proportional to 7fco, implying a small vari- 
ation oi A/B for small 7. 

Furthermore, the form p4(l is expected to be valid in 
the regions where the variation A/B is small. By looking 
at fig. El with the axes interchanged, one recognizes that 
the variation of A/B is strong near the gap edges, but 
weak otherwise. 

Ignoring this potential shortcoming in the ansatz (|14(l , 
and neglecting terms involving S"{w), the determinental 
consistency equation (|15() for A and B reduces to a 4th 
order algebraic equation for S'^w). This equation can be 
solved analytically, giving a very long and complicated 
expression S'l for this approximation to S'{w). 

In fig. we compare this solution S[ to a further sim- 
plified approach: after switching back to the old z variable, 
we take 7 — > to obtain S2, this corresponds exactly to 




Fig. 5. Comparison of the two approximations made: we plot 
the relative error (l^i — S'2|)/|5'i| between the two approxima- 
tions 5*1 and 5*2 over the z-range. The two lines correspond to 
the two solutions, see the two branches in fig. |3| In the case of 
the branch with the forbidden zone (contirmous line), the ap- 
proximation 1161 for obtaining 52 is valid a long distance away 
from the points 21 and 22, but breaks down in their vicinity. 

the de Vries wavevector function derived from eq. HlOfl and 
used in fig.|31 S2 stands for the local wave vector k{z) of a 
plane wave with frequency w = fepc in a cholesteric with a 
inverse pitch q{z), which we consider as locally constant - 
thus taking 7 = here means that we ignore derivatives of 
the phase (/>, in other words an adiabatic approximation: 

5^-±(iA:2(ei+e2) + g2 (le) 

±^ko^{ei - e2)2fc^ + Sq^ei + £3)) 

where q — q{z) again is given by eq. H12|l . 

We note that, for the branch with the forbidden zone, 
the approximation which takes us from S'l to S2 is valid 
far away from the points zi and Z2, where the phase varia- 
tion is zero on the transition from being real to becoming 
imaginary or vice versa. However, close to these points, 
approximating S'l with S2 clearly breaks down, as fig. |S1 
shows. In any event at the band edges where S'{z) the 
WKB-like approach breaks down for the usual reasons. 

4.4 Comparison of approximate schemes with 
numerical solutions 

So far, we have only compared the validity of two consecu- 
tive approximations S'l and S2 with each other. However, 
the true test is to compare the approximation to the full 
numerical solution. Figure compares the numerical solu- 
tion with the approximate solution obtained by using 
For this purpose, we choose the approximate solution ^2 
corresponding to the branch (2) in fig. |3 integrate it nu- 
merically and obtain the resulting wave by eq. (|14|l . The 
numerical solution is obtained by first evaluating eq. (|14|l 
at an anchor point zq giving the four necessary initial con- 
ditions: the two field components and i?^, and their 
gradients E'^ and E'^^. In fig. El (a), we chose zq = 0, and 
in figs. El (b)-(d) Zq = 70. We observe that the agreement 
between the two solutions breaks down on approaching 
zi and Z2 ■ This has a plausible explanation: we have seen 
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Fig. 6. Comparison of the full numerical solution (continuous thin line) and the approximate solution (overlaid dashed line): 
the approximate solution is based on an numerical integration of the local de Vries expression 1161 and eq. 1141 . In fig. (a), the 
numerical solution is forced to have the same initial conditions at the point zo — 0, whereas in figs, (b)-(d) zq — 70 (circled). 
Plots (c) and (d) demonstrate the stability of the approximate solution for values of positions z distant from zq. Note that the 
2;-range is different in each of the diagrams (b)-(d) indicating faster oscillations of the field {E^,En) for large z in the same way 
as q{z) increases according to eq. I|12^ . 



in fig. that the approximation S2 breaks down near the 
points zi and Z2 as in WKB. We therefore expect the cor- 
responding solution of the wave equation that is generated 
from 5*2 to fail dramatically at the same point. 

The discrepancy at zi and Z2 between the approximate 
and numerical solution can also be explained from an- 
other viewpoint: between zi and Z2, we ideally expect two 
fundamental solutions, the one increasing and the other 
decreasing in an exponential-like fashion. As the proce- 
dure described in the last paragraph introduces small er- 
rors in the initial conditions at zq, the numerically com- 
puted solution does carry components of the two funda- 
mental solutions. As z enters the region between zi and Z2, 
one solution decays rapidly, but the other (initially only 
marginally present) increases quickly, becoming the domi- 
nant contribution in the region [zi , Z2] ■ Away from the gap 
the agreement is very good - over large distances above 
the gap (indicating great stability) and even backwards to 
the point z = —20 where locally the pitch is infinite. 



4.5 Transmission coefficient 

One can estimate the dependence of the transmission co- 
efficient on the principal variables qo, 7, ko, ei and €2- 



Our system is analogous to the simple quantum me- 
chanical barrier problem for a free particle (see for exam- 
ple |H] ; here we have a vector field E with two components 
rather than one): outside a region [0, 1], propagation is de- 
scribed by the Schrodinger equation of a free particle with 
an oscillating wave function and allowing for an incoming 
and an outgoing wave. 

However, in the region [0, 1], the particle moves across 
a potential V. If the potential V is constant and greater 
than the particle's energy E, the solutions to the wave 
equation are either exponentially decaying or increasing. 
The reflection or transmission coefficients can easily be 
calculated exactly. For a general potential barrier V{z), 
the transmission coefficient can be estimated jQ: 

T«exp|-2i /" dz^y2m{V{z) - E)\ . 



We estimate the transmission coeflicient T as the square 
of the ratio of the fields or E,j at zi and Z2: 
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Fig. 7. Quality of the approximation l|18|l to model the func- 
tion 52 in the evanescent region between zi and Z2: the thin 
curve corresponds to the imaginary part of S'2, whilst the over- 
laid dashed curve is the approximation according to eq. 11811 . 
The inset shows the relative error. 



using eq. ((T^ for the scaling of 21 — Z2 . 

To arrive at this estimate (|19f) of T, we used the ap- 
proximation £"2(2), which breaks down near the edges of 
the stop gap at zi and 22, as fig. |S1 shows. However, as 
one can see from fig. [7| the dominant behaviour of T is 
determined by the large imaginary values of S'2 ■ These oc- 
cur deep inside the stop gap where the approximations are 
valid. 

Our expectation at eq. (|13|l has been confirmed: the 
wider the spatial extent of the stop gap the smaller the 
penetration and hence the lower the transmission. Small 
7, that is slow spatial variation of the band structure, gives 
lower transmission (however providing that the sample is 
thicker than the spatial extent of the gap!). One requires 
7 ^ {n2 — ^1)^(^0/(70)^ for effective reflection at the ef- 
fective frequency k^. 



where the latter form arises from our assumption (|14l) . 
The barrier form of T now yields 



exp(2z[£(z2)-5(zi)]) 
s exp [ 2i dz 5*2(2) 



(17) 



where we have implicitly assumed that the integral of 
5*2(2) over z between zi and Z2 is purely imaginary, yield- 
ing a real number in the final exponential of eq. H17|l . In 
the forbidden zone [zi , 22] , there are two solutions for S'2 
(see. eq. (jl6|l ^: both are purely imaginary, but of opposite 
signs. Therefore, we have to choose the correct one to be 
inserted into eq. H17|l . corresponding to the decaying wave. 

The function 5*2(2) in (|16|l has divergent slopes at the 
boundary points Z2 and zi and assumes a maximum in 
between with 



S' 



£2 - ei 



62 



Therefore one might approximate 52(2) by an ellipse be- 
tween the points zi and Z2 with a maximum value of 



^2 max" 



'^2app(^) 



(18) 



£1 + £2 



[q - V^iko) (%/e2fco - q) 



where q = q^ + q^^z. Figure[7|shows the £2(2) and overlaid 
the approximation 5*2 app of eq. H18|l . The inset shows the 
relative error (5*2 - 5'2app)/5'2. 

Finally, this expression can be integrated analytically 
to give an estimate for the integral in eq. (|17|l . which cor- 
responds to the area under the curve in fig. |7| We obtain 



T exp 



exp 



TT kl 



\/24<?g7 Vei + £2 

21 \2 \ 

%1 (^1 - \ 



V24 VeT 



(2 



(19) 



5 Conclusions 

The numerical solution of the field amplitudes for light 
in a spatially varying cholesteric (a photonic solid) has 
been presented. Fields decay in the region where, for the 
relevant frequency of light, a uniform cholesteric would 
have its stop gap. Transmission is low for sufficiently slow 
variation of the pitch. Thus the simple idea that each com- 
ponent of a broad spectrum of light is reflected where its 
Bragg (strictly its de Vries) condition is satisfied, is par- 
tially borne out. 

Approximate field amplitudes are derived by taking 
this local de Vries condition literally. They are reasonable 
approximations far from the gap, but less good at the gap, 
for reasons we discuss. 

Slower spatial variation, 7 of pitch gives much lower 
transmission, in fact exponentially lower, varying as e~'^/'^ . 
This leads to better utility as broad band reflectors in de- 
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